How Far Can the SO (10) Two Higgs Model 



Describe the Observed Neutrino Masses and Mixings ? 
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Can the SO(10) model with one 10 and one 126 Higgs scalars give the observed masses and 
mixings of quarks and leptons without any other additional Higgs scalars? Recently, at least, for 
quarks and charged leptons, it has been demonstrated that it is possible. However, for the neutrinos, 
it is usually said that parameters which are determined from the quark and charged lepton masses 
cannot give the observed large neutrino mixings. This problem is systematically investigated, and 
it is concluded that the present data cannot exclude SO(10) model with two Higgs scalars although 
it cannot give the best fit values of the data. 



PACS number(s): 14.60. Pq, 12.15.Ff, 12.10.-g 



I. INTRODUCTION 

SO (10) GUT model seems to us the most attractive 
model when we take the unification of the quarks and 
leptons into consideration. However, in order to repro- 
duce the observed quark and lepton masses and mixings, 
usually, a lot of Higgs scalars are brought into the model. 
So it is the very crucial problem to know the minimum 
number of the Higgs scalars which can give the observed 
fermion mass spectra and mixings. A model with one 
Higgs scalar is obviously ruled out for the description of 
the realistic quark and lepton mass spectra. Two Higgs 
models were initially discussed by Mohapatra et.al [Q. 

In the previous paper we discussed 2 Higgs scalars, 
{10 and 126} case and {10 and 120} case, and showed 
that they reproduce quark-lepton mass matrices unlike 
the conventional results ||. One of new points of our 
approach is that we adopt general forms of Yukawa cou- 
plings allowable in the SO(10) framework. However, we 
did not argue there about the neutrino mass matrix since 
it may incorporate additional assumptions like the seesaw 
mechanism etc. 

One of the merits of the SO(10) model is that it in- 
cludes a right-handed Majorana neutrinos in the fun- 
damental representation and naturally leads to the see- 
saw mechanism. Also some papers claim that two Higgs 
model ({10 and 126+126}) does not reproduce the large 
mixing angle of the atmospheric neutrino deficit Jij . So in 
this paper we apply our method developed in the previous 
paper to the neutrino mass matrix, fitting the other pa- 
rameters of the quark-lepton mass matrices. Our model 
has the two Higgs scalars {10 and 126} both of which are 
symmetric with respect to the family index. Therefore 
mass matrices are symmetric whose entries arc complex 
valued. We do not adopt another choice {10 and 120}. 
For it does not involve the mass term of the right-handed 



Majorana neutrinos which are the ingredients of the see- 
saw mechanism. 

We begin with the short review of our previous work 
0. In the case where two Higgs scalars, 4>iq and 0126, 
are incorporated in the SO(10) model, the mass matrices 
of quarks and charged leptons have the following forms 



M, 



c M + dMi, M d = M + Mi, M e = M 



3Mi. 
(1.1) 



Here Mo and Mi are the mass matrices generated by the 
Higgs scalars <fiia and </>i26, respectively. Also cq and C\ 
are the ratios of VEV's, 



(1.2) 



and <j> u and <fi d denote Higgs scalar components which 
couple with up- and down- q uar ks, respectively. Elimi- 
nating Mq and Mi from Eq.([hl]), we obtain 



where 



Cd 



M e = c d M d + c u M u , 



3c + ci 



c - ci 



c - ci 



(1.3) 



(1.4) 



Since M„, M d , and M e are complex symmetric matrices, 
they are diagonalized by unitary matrices U u , U d , and 
U e , respectively, as 



UlM u U u = D u , UjM d U d = D d , UjM e U e = D e 



(1.5) 



where D u , D d , and D e are diagonal matrices given by 
D u = diag(m„,m c ,TO t ) , D d = di&g(m d ,m s ,m b ) , 



D e = diag(m e , m M , m T ) 



(1.6) 
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Since the Cabibbo-Kobayashi-Maskawa (CKM) matrix the relation ([L3) is re-written as follows: 
V q is given by 



(1.7) 



{UtU u ) T D e (UtU u ) = c d V q D d V q T + c u D u . (1.8) 
Therefore, we obtain the independent three equations: 



detD e Dl = \c d \ 6 det 



TrD e Dt = \c d \ 2 Tr[(V q D d VT + KD u )(V q D d V^ + kD u )^ 
Tr(D e Dl) 2 = \c d \ 4 Tr\((V q D d VT + KD u ){V q D d V^ + kD^) 2 



(V q D d V q T + K D u )(V q D d V q T + kD u )^ 



(1.9) 
(1.10) 
(1.11) 



where k = c u /c d . By eliminating the parameter c dl we 
have two equations for the parameter k: 
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2(m 2 m 2 + m^m 2 . 



%m\) (|1.9|) 2 - (|1.1C|) 



(1.12) 
(1.13) 



whe re (|l.9[ ) 3 , for instance, means the right-hand side of 
Eq.(1.9) to the third power. Let us de note th e par ameter 
values of k evaluated from Eqs.(LlS) and ( 1.13 ) as ka 
and Kb, respectively. If ka and kb coincide with each 
other, then we have a possibility that the SO(10) GUT 
model can reproduce the observed quark and charged lep- 
ton mass spectra. If ka and kb do not so, the SO(10) 
model with one 10 and one 126 Higgs scalars is ruled out, 
and we must bring more Higgs scalars into the model. 

Note that Eqs. (1.9)- (1.11) can constrain only the ab- 
solute value of c d = |cd|e l<T . The argument of the param- 
eter c d can be determined by taking neutrino sector into 
consideration. In the previous paper [jO, we have found 
that only for the signs of the masses 

(m t , m c , m u ;m b , m s ,m d ; m T ,m^, m e ) 

= (+,-,+;+,-,-; +,±,±) (a), (1.14) 

and 

= (+,--;+,-, -;+,±,±) (b), (1.15) 

there are solutions which gives ka = kb, and the corre- 
sponding parameter values (|cd|,«) are 

(\c d \, k) = (3.15698, -0.019296e 2 64172Oi ), (1.16) 



(3.03577, -0.019398e 2 " 57 ° Oi ) for (a), (1.17) 



and 



= (3.13307, -0.019314e 2 n464Oi ), (1.18) 
(3.00558, -0.019420e 310014Oi ) for (b) (1.19) 

and m s = 76.3 [MeV] for input 6> 23 = 0.04 20 [rad] and 
8 = 60° at /j. = mz (mz is the neutral weak boson mass) . 
For the relation between the values at fj, = mz and those 



at (j, = Ax (Ax is a unification scale), see Ref. 0. The 
purpose of the present paper is to investigate whether 
these solutions can give reasonable values for observed 
neutrino masses and mixings or not. 



II. THE NUMBER OF PARAMETERS IN THE 
SO(10) MODEL WITH TWO HIGGS SCALARS 

As we have discussed in the previous section, among 
four freedoms of complex {co, ci} or {c dl k}, we have 
been able to fix the three of them, k and \c d \. This is not 
accidental. Let us discuss the situation in detail in the 
SO(10) two Higgs model. 

In the previous paper, by using the relation ([ll]), we 
have investigated whether there is a set of parameters 
which can give the 13 observable quantities D e , D Ul D d , 
and V q or not. We can rewrite Eq.(1.8) as 



where 



A T e D e A e = c d (V q D d V q T + kD u ), 



A e = U\U U , 



Cd e 



(2.1) 



(2.2) 



(2.3) 



The quantities D e , D u , D d , and V q are inputs, and 
the quantities \c d \, k, and A e are the parameters which 
should be fixed from those observed quantities. In gen- 
eral, an n x n unitary matrix for n generations has n 2 
parameters. Therefore, the number of the parameters is 

AT(pmt) = N(A e ) + N(c d ) + N(k) =n 2 + 2 + 2. (2.4) 

On the other hand, the number of equations is 

N(eqs) = n(n+ 1), (2.5) 



because Eq.(2.1) is symmetric. Therefore, the number of 
the unfixed parameters is given by 



AT frcc = iV(pmt) - TV(eqs) = 4 - n = 1, 



(2.6) 
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for n = 3, i.e., the 13 observed quantities fix the param- 
eters \cd\, k, and A e , but 1 parameter a remains as an 
unknown parameter. 

In the present paper, we will try to predict neutrino 
masses 

D v = U^M U U V , (2.7) 
and mixing matrix 

V t = U*U* U , (2.8) 

by using the observed quantities D e , D u , Dd, and V q a nd 
the parameter values \cd\, k, and A e fixed by Eq.(2.1). 

SO(10) GUT asserts that the Dirac neutrino mass ma- 
trix Mfl is given by the form 



N{A U ) = N(Vi) = n 2 



(2.21) 



M D = c M - 3ciMi, 



(2.9) 



and Majorana mass matrices of the left-handed and right- 
handed neutrinos, Ml and Mr, are proportional to the 
matrix M\\ 

M L = c L M 1 , M R = c R M x , (2.10) 

where Mq and Mi are related to the quark and charged 
lcpton mass matrices M u , Md, and M e as follows: 

3M d + M e 



Mn = 



Mi = 



M d -M e 



(2.11) 
(2.12) 



Then the neutrino mass matrix derived form the seesaw 
mechanism becomes 

M v = M L - M D M^Ml 
= c L M 1 

-c^{cqM q - 3c 1 M 1 )Mi 1 (c M - 3c 1 M 1 ) T . (2.13) 

In the present paper we adopt cl = 0. Also we may ig- 
nore the phase of cr which does not affect the observed 
values. Therefore, we can rewrite Eq.(2.13) as 



\c R \AlD v A v = MdM^MI, 



similarly to Eq.(2.1), where 



M D = c Q M Q - 
Mo = i(3M d 

Ml = i(M d - 



3ciM 1; 
+ M e ), 

■M e ), 



with 



M d = UlM d U u = V q D d V q T , 
M e = U^M e U u = A T e D e A e 



Cd(V q D d V* 



kD v 



(2.14) 

(2.15) 
(2.16) 

(2.17) 

(2.18) 
(2.19) 



Differently from the previous work, the quantities D u and 
Vg are unknown parameters at the present stage. Since 

Vt = A* e Al, (2.20) 



and A e is fixed fro m Eq. (2.1), the number of the unknown 
parameters in Eq.(2.20) is 



Of course, the unknown parameters in A v contain the n 
unphysical parameters which cannot be determined be- 
cause of the rephasing in the fields e^. Therefore, the 
number of the unknown parameters is 

TV(pmt) = N(D V ) + N{A U ) + N(\c R \) + N(a) 

= n + n 2 + 1 + 1 = n 2 + n + 2 (2.22) 

and from the number of equations TV(eqs) = n(n + 1) in 
Eq. fl2.14 ), we obtain the number of the unfixed parame- 
ters as 



TV(pmt) - TV(eqs) 
(n 2 + n + 2) — n(n - 



1) 



(2.23) 



This means that we can predict neutrino masses and mix- 
ing completely if we give the two values |cr and a. The 
numerical predictions will be investigated in the next sec- 
tion. 



III. NUMERICAL RESULTS 

Here we discuss the numerical results of the neutrino 
mass spectrum and neutrino mass matrix. For exam- 
ple, we use the set in Eq.(1.18). Even if the other sets 
are used, our results are scarcely changed. The allowed 
values of neutrino mass square differences and lepton fla- 
vor mixing angles depict complicated tracks with moving 
a = arged (Fig. 2). This figure shows a general ten- 
dency that the lepton flavor mixing angles #12 and #23 
get larger as a approaches to 37r/2. For an illustration 
we take a = 1497r/100, then these values become 



A? 



'12 



l 13 



0.15, 



sin 2 (26»i 2 ) = 0.76, 
sin 2 (26» 13 ) = 0.16. 



0.85, 



Am| 3 
Am 2 3 



sin 2 (20 23 ) = 0.75, 



(3.1) 



There still remain a little bit discrepancies between our 
results and experiments. However our results are much 
improved in comparison with those by Babu-Mohapatra 
Q in which they obtained sin 6>i2 = — 0.3, sin $13 = 0.05, 
and sin #23 = 0.12 — 0.16. The purpose of the present pa- 
per is to study the general tendency of the fittings and 
not to pursuit the precise data fitting, for the data them- 
selves are not affirmative and we have theoretical ambi- 
guities not incorporated in the present data fitting like 
the renormalization gr oup effect. 
In the choice of Eq.(3.1), we have 



|c d | = 3.16 
co = 



= 54.84e 



-20.24°i 



Cl = = 70.54e +41 - 90 °\ 



(3.2) 
(3.3) 

(3.4) 



In this case, Eq.( |2.11 ) - E q.(2.15 ) are re- written in the 
basis of M u = D u (sec Eq.(l.Sj]as 



3 



3V q D d V q T + c d {nD u + V q D d V q T ) 
Mo = 

/ -0.00405e~ 57 - 29 ° l -0.00753e~ 56 - 24Oi -0.00533e+ 65 - 46 ° 4 
= 2.1646 x io 3 e + 10 - 48 ° l -0.00753e~ 56 - 24Oi -0.02986e~ 51 - 59 ° J +0.06358e" 57 - 64 ° l ) [MeV], (3.5) 

V -0.00533e +65 46 ° l +0.06358e- 5764 ° 4 +1.00000 

V q D d V q T - c d {nD u + V q D d V q ) 
Mi = 

/ -0.00715e+ 95 23 ° l -0.01333e +96 54Oi +0.00944e +38 23Oi 
= 9.5127 x 10 2 e- 24 44 ° l -0.01333e+ 96 - 54 ° 4 -0.04878e+ 90 - 73 ° l +0.11247e+ 9513 °* ) [MeV] , (3.6) 

V +0.00944e+ 38 - 23 ° 4 +0.11247e +9513Oi +1.00000 

\cr\M v = (c M - SciAfiJMf^coMo - 3 Cl M!) T 

+0.1163e +26 - 89 ° l +0.2165e+ 28 06Oi -0.1536e- 3a53 ° l 
-4.6628 x 10 6 e- 5217 ° l ( +0.2165e+ 28 06 ° l +0.8193e+ 28 00 °* -1.9276e +29 - 52 ° l | [MeV]. (3.7) 

-0.1536e- 30 - 53 ° 4 -1.9276e +29 - 52 °* +1.0000 

Let us choose the free parameter \cr\ so as to result in small neutrino masses, for example when |c^| = 3.2 x 10 14 , we 
have Am 2 ^ = 1.5 x 10~ 3 [eV 2 ]. 

Here there arises a question what makes the two flavor mixing angles large. We need to investigate the mixing 
matrices U e and U v which diagonalize M e and M v , respectively. Those are obtained as 

-0.863 +0.504e+ 9 46 ° 4 -0.022e +56 - 66 ° 4 

U e = ( -0.493e" 9 82 ° l +0.834 -0.248e +16 - 63 ° 4 ) , (3.8) 

-0.110e- 21 - 40 ° l +0.223e- 1810 ° 4 +0.969 

+0.992 -0.092e- 15 - 94 ° l -0.088e+ 12 86 ° l 

U v = ( +0.049e+ 76 - 86 ° l +0.724 -0.688 e - 16 08 ° 4 ) . (3.9) 

-0.117e+ 9 80 ° l +0.683e +16 74 ° l +0.721 

Here, |t/ e ii|, C4i2|, \U e 2i\, \U e 22\ > 0.5 for the charged lepton mass matrix and f/^ 2 2 1 , \U U 2z\, \UuS2\, \U v a&\ > 0.7 for 
the neutrino mass matrix. Therefore the components of the lepton flavor mixing matrix become |Vqi|, |Vji2|, |Vz2i|, 

|^22|, |Vj 23 |, \Vl 32 \, 1^33 1 £0.5: 

-0.844e+ 210 ° l -0.494e" 9 95 ° 4 +0.206e +23 - 61 ° l 

Vi = | +0.527e+ 3 26 ° l +0.696e" 8 84 ° 4 -0.488e+ 24 - 970 * | . (3.10) 

-0.098e- 15 - 78 ° l +0.521e- 27 - 43 ° l +0.848e+ 6 - 32 ° 4 

The mixing angle 623 becomes larger, while the mixing angle $12 smaller, if we take the smaller value of \mt\ or \m d \, 
or larger |m c |, or \mb\, \m s \ than their center values. 

As a simple example, the shift of \m d \ and \m s \ causes the change of mixing angles and neutrino mass square 
differences as depicted in Fig. 2. Fig. 2 shows that the 023 and $13 can appr oach the 99%C.L. of SK §+CHOOZ @ 
but 6»i2 and Am 2 2 are out of the range of 99% - 99.9%C.L. of SOLAR ff+CHOOZ. 



4 




FIG. 1. The relation between our results and the two-flavor oscillation analysis [6] when a is moved, (a) The circles and 
triangles indicate the values of Am 2 3/Am 2 3 and Ami 2 /Ami3 at every 7r/2 of a. (b) The circles, triangles, and stars indicate 
the values of sin 2 2023, sin 2 2#i2, and sin 2 2#i3 at every 7r/2 of a, respectively, (c) The circles, triangles, and stars indicate the 
values of (Am 2 3 , sin 2 2^23 ), (Amf 2 , sin 2 2#i2), and (Am 2 2 , sin 2 2^13) at every ty/2 of a. Here we have set Am^ = 1.5 x 10~ 3 [eV 2 ] 
in every case. 
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12 

FIG. 2. The relation between 3^ Oscillation analyses by G.L.Fogli et al. [10] and by us for Am| 3 = 1.5 x 10~ 3 [eV 2 ]. (a 



For SK+CHOOZ. (b) For SOLAR+CHOOZ. The circles indicate our solutions for Eq.(^lj). The solid line through them is the 
track as m<j is varied. It goes from the experimental limit that \m,d\ moves over the range, 4. 03 - 5.29 MeV_[5]. At that time, 
\m s \ simultaneously changes over the range, 76.3 - 76.2 MeV so as to satisfy the relations ( fl.l2| ) and (1.13). If we take the 
smaller \m,d\ with the fixed <r, the solution in (a) moves rightward and the solution in (b) does left-upward (Table I (i)). Since 
the minimum \rrid\ for (b) gives bad fitting, we have changed o from 1497r/100 to 1467r/100 , which is denoted by star (Table I 
(ii)). Thus our result approaches the 99%C.L. of SK+CHOOZ and 99.9%C.L. of SOLAR+CHOOZ. 



IV. DISCUSSION 



Since there are only two basic matrix Mo a nd M i in 
this model, the number of parameters in Eq.( [2.l| ) and 
(f2T4l) is 



D u , D d , D e , D w 3x4 = 12 
c d , |cr|, k 2+1+2 = 5 

V q , Ae, A v 4+9+9 = 22 
sum. 39 



(4.1) 



and the number of equations is -/V(eqs) = 12 x 2 = 24. 
Therefore the number of free parameters is iV(pmt) — 
jV(eqs) = 39 — 24 = 15. On the other hand, the number 
of the physical parameters which can be determined by 
experiments is 



m u , m c , m t , 




3 


ma, m s , m>b, 




3 


CKM: 9 12 , 6 23 , 




4 


m e , rrift, m T , 




3 


m Vs , m Vtl , m VT , 




3 


MNS: 9 12 , 9 23 , 


013) P> P 


6 


sum. 




22 



(4.2) 



where (3 and p are Majorana phases in the Maki- 
Nakagawa-Sakata (MNS) matrix because of no rephasing 
in the neutrino fields vi,. To sum up the matter, we dis- 
cuss the consistency test about 22 physical parameters 
by using only 15 free parameters. The consistency test 
in the quark sector is good, as shown in our previous pa- 
per. In the lepton sector, the test is not so bad when we 



G 



adopt the MSW large mixing angle solution of solar neu- 
trino deficit, and this model favors the normal hierarchy 
of neutrino mass spectrum. 

Also we can predict the yet unobserved values such as 
the averaged neutrino masses (m) a p and Jarlskog param- 
eter in the lepton part. The averaged neutrino masses 
appear in the reactions where the Majorana neutrinos 
propagate in the intermediate states. They are 



(m v )af. 



3 

E 



U aj U, 



3 U Pj> 



(4.3) 



where a and f3 are (e, /x, rj 
noless double beta decay fLl 



They correspond to neutri- 
for a = (3 = e, fi — e conver- 



sion (/j,~ + (A, Z) 
K decay (K ~~ — > 



f {A, Z -2) tax a 
for a = (3 = 




Fig. 3 we have depi cted a dependence of (m v ) a p / \/K 



In the case of Eq. 
(m) a p 



°23- 



1|), these values become as follows. 
0.032 



'23 



0.065 0.30 
0.096 0.59 
0.67 



(4.4) 
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For instance, if we input Am| 3 = 1.5 x 10~ 3 [eV 2 ], (m) ee 
becomes 0.0012[eV]. This is smaller than the experimen- 
tal value of the next generation experiments such as GE- 
NIUS pi, CUORE Q, and MOON @. Jarlskog pa- 
rameter^! appears in three generations 

n / \ i->/ \ 7 A£ , 2iA£ , 32 Ai?3i 

P{v e -» z^J - Piv^ -> u e ) 



J 



xsm 



sin 



fAE^L 



sin 



( AE$[L 



with 

Here we have adopted the notation 



AEjk = Ej — Ek — 



Ar 



2£ 



M 



(4.5) 



(4.6) 



(4.7) 



with 



f/diag(£i, E 2 , £3)^ + diag(a, 0, 0) 



[/ M diag«, 



J5. 



M 



){U M ) 



M\-l 



(4.8) 
ID, 



^2 ,-^3 

The cr dependence of J is depicted in Fig. 4. For Eq. 
it takes 

J ~ 0.00015. (4.9) 

However, it needs careful consideration that J drastically 
changes at a ~ 37r/2. (m) a p and J in the cases of Table 

I (i) and (ii) discussed in Fig. 2 are also listed in Table 

II (i) and (ii). In this paper we have discussed how far 
the SO(10) two Higgs scalar model describes the quark- 
lepton masses and mixing parameters. We can conclude 
that this model cannot be rejected within the existing 
data. It should be remarked that the whole parameters 
can be decided from the existing data in principle. 
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FIG. 3. The relations between the averaged neu- 
trino masses of lepton number violation process and a 
The hexagons, white circles, boxes, triang les, bl ack cir- 
cles, and stars indicate the values of (m) e 

y/ Am 23> ( m )"/\/ A ™23. Ww/\/ Am 2 



/-y/Amfa, (m) eM / 
, (m) MT / v /Am| 3 , 



and {m) TT / \J Am^ at every 7r/2 of cr, respectively. 
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FIG. 4. The relation between Jarlskog parameter J and 
a. The circles indicate the values of J at every 7r/2 of a. 
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(ii) 



\m d \ =4.03[MeV] 
(Amf 2 )/(Am? 3 ) = 



2 (20i 



\m s 
0.43 

,2 



76.3[MeV], o = 149tt/100, 
(An4)/(Am? 3 ) = 0.57, 
) = 0.17 



0.52, sin 2 {2023) = 0.91, sin 2 (20 



\m d \ =4.03[MeV 
(Amf 2 )/(Am 13 ) = 0.20, 
sin 2 (20i 2 ) = 0.54, sin 2 (20 ; 



m s \ = 76.3[MeV], o = 146tt/100, 
(Ami 3 )/(Am 2 3 ) = 0.80, 
= 0.88, sin 2 (20i 3 ) = 0.20 



TABLE I. Our solution (the second and third lines) from 
the input parameters (the first line). The result of (i) is ob- 
tained when we move \m d \ from 4.69[MeV] to 4.03[MeV]. (ii) 
is the result when we move \nid\ as (i) and, furthermore, a 
from 1497r/100 to 1467r/100. These data fitting corresponds 
to Fig.2. 



(0 


H«/V Am 23 


= 0.039, 


(m) eM / v /Am 2 3 


= 0.086, 






= 0.43 , 




= 0.19 , 






= 0.96 , 


W"/\/H 


= 0.52 , 




J = 0.0091 








(ii) 


(m} ee /-y/Am| 3 


= 0.028, 


(m) e n/y/Am% 3 


= 0.064, 




(m) eT /-v/Am:jg 


= 0.32 , 




= 0.10 , 






= 0.71 , 


(m) TT / v /Am| 3 


= 0.54 , 




J = -0.014 









TABLE II. The values of averaged neutrino masses and 
the Jarlskog parameter for the case (i) and (ii) in Table I. 
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